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The dynamical generation of entanglement in closed bipartite systems is investigated in the semi- 
classical regime. We consider a model of two particles, initially prepared in a product of coherent 
states, evolving in time according to a generic Hamiltonian, and derive a formula for the linear 
entropy of the reduced density matrix using the semiclassical propagator in the coherent-state rep- 
resentation. The formula is explicitly written in terms of quantities that define the stability of 
classical trajectories of the underlying classical system. The formalism is then applied to the prob- 
lem of two nonlinearly coupled harmonic oscillators and the result is shown to be in remarkable 
agreement with the exact quantum measure of entanglement in the short-time regime. An im- 
portant byproduct of our approach is a unified semiclassical formula which contemplates both the 
coherent-state propagator and its complex conjugate. 
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I. INTRODUCTION 

Entanglement is one of the most formidable effects 
of the quantum world. Its puzzling nature, which in- 
trigued the scientific community for a long time, is now 
being used to accomplish tasks such as quantum infor- 
mation processing, quantum computation, teleportation 
and quantum cryptography [H, Also, its importance 
has been recognized in the context of several foundational 
issues underlying the quantum theory, from the explana- 
tion of the quantum-classical transition — and its im- 
plications to the measurement problem — to the under- 
standing of the nonlocal aspects permeating the EPR 
debate g-i]. 

Entanglement is widely believed to be a purely quan- 
tum effect with no classical analogue. Despite this com- 
mon belief, several results have been reported associating 
the entanglement dynamics with classical quantities. For 
instance, in Refs. p-Q it is shown that the entanglement 
dynamics can be approximately simulated in the short- 
time regime by the Liouvillian formalism. In particular, 
for some specific couplings, the Liouvillian entropy has 
been shown to reproduce exactly the entropic measure 
of entanglement for all values of time @. In addition, 
in Ref. [8[ the authors have analytically shown that the 
short-time dynamics of entanglement does not depend 
on h for a large class of Hamiltonian systems. Finally, 
some authors investigated entanglement in the semiclas- 
sical regime by means of time-dependent perturbation 
theory QUEI. 

The scenario delineated by these works points to a sit- 
uation in which a statistical theory based on classical 
trajectories is able to predict the dynamics of a quan- 
tity meant to be exclusively quantum [l2| ■ This observa- 
tion leads us to suspect that the entanglement dynamics 
is initially promoted by mechanisms with well defined 
classical analogues. Finding out these mechanisms is the 
main motivation of this paper. We follow, however, a pro- 
gram that is substantially different from the works quoted 



above, as it is based on semiclassical methods instead of 
classical statistical theories. Specifically, we propose to 
derive a semiclassical measure of entanglement in terms 
of the semiclassical propagator in the coherent-state rep- 
resentation [l^ - tl7j . 

Recently, a similar calculation has been carried out [Tel 
Il9| which differs from ours in some important aspects. 
First, the approach adopted there was based on the Van- 
Vleck semiclassical propagator [2(| , 



E 

traj. 



1 d 2 S(q 2 , qi ,T) 



2-kH dq^dqi 



,j-S(q 2 ,qi,T) 



This is a semiclassical formula for the one-dimensional 
quantum propagator {q2\e~ lHT / h \qi) in the position rep- 
resentation. This formula depends only on the classical 
trajectories of an underlying classical dynamics connect- 
ing the initial position q\ to the final position q2, dur- 
ing the time interval T. The function S(q2, qi,T) is the 
classical action of the trajectory and the sum runs over 
all trajectories satisfying the boundary conditions. Our 
approach, on the other hand, is based on a semiclassi- 
cal propagator formulated in the coherent-state repre- 
sentation, which has the advantage of offering a straight- 
forward extension to systems with spin degrees of free- 
dom [HH23. 

A second important difference relics on the fact that 
our approach does not employ any averaging over the 
initial conditions in phase space. Even though in 
Refs. [l8|, this statistical procedure is claimed to be 
nonrestrictive we believe it is not mandatory from a phys- 
ical point of view. The only approximations used here are 
those usually associated to the method of the stationary 
phase. 

Finally, we observe that, contrary to the Van-Vleck 
propagator, the coherent-state propagator is generally 
determined by complex trajectories and actions. This in- 
troduces an additional technical difficulty, namely, that 
the complex conjugate of the semiclassical propagator 
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does not have a straightforward interpretation. Actually, 
this turns out to be an interesting mathematical issue 
to be understood in the context of general applications 
of the coherent-state propagator. In this paper we for- 
mulate and address this problem as a preliminary step 
towards the derivation of a semiclassical formula for the 
entanglement. 

This article is organized as follows. In Sec. HU we 
present the main ingredients of the original formula of 
the semiclassical propagator in coherent states and then 
extend it to contemplate also the complex conjugate of 
the propagator in a unified formalism. We then proceed 
to calculate the semiclassical purity in Sec. IIII| where a 
further extension of the semiclassical formula of the prop- 
agator is required to accommodate bipartite systems. In 
Sec. |lV]we present a study of case; the formalism is ap- 
plied to the problem of two nonlinearly interacting oscil- 
lators. Finally, in Sec. [V] we summarize and conclude the 
work. 



II. SEMICLASSICAL PROPAGATOR IN THE 
COHERENT-STATE REPRESENTATION 

The aim of this section is twofold. Firstly, we briefly 
review some of the main aspects of the semiclassical for- 
mula of the coherent-state propagator. For a subsid iary 
literature on this representation we refer to Refs. |26|-[28|. 
Secondly, we show how to extend the formalism so as 
to semiclassically approach both the propagator and its 
complex conjugate in a unified mathematical structure. 
In this sense, our approach intends to offer a generaliza- 
tion of the formula derived in Ref . [l?} . 

We start, as a motivating question, with the gen- 
eral problem of calculating the expectation value of an 
arbitrary operator A via semiclassical propagator. If 
the initial state of the system is the coherent state 
| zo) , then at the instant T the mean value (A)t = 
{zo\e lHT / h Ae~' lHT / h \z ) can be written in terms of the 
propagator as 

f d 2 z 2 d 2 Zl , 
(A) T = / — -K(z ,z 1 ,-T)A{z 1 ,z 2 )K(z2,z ,T), 



where A(z\,Z2) — {z\\A\z2) and 

K(z ,zx,-T) = (z \e ii ' T / H \z 1 } = K*(z 1 ,z ,T). 

The complex conjugate of the propagator K* is going to 
be present whenever measurable quantities are regarded. 
However, to the best of our knowledge, there is no pre- 
scription on how to obtain the semiclassical version of this 
object in the coherent-state representation. But should 
not we simply take the complex conjugate of the semiclas- 
sical propagator? We opt here by a more careful strat- 
egy that preserves both the interpretation of the critical 
trajectories and the rigor of the original derivation in 
Ref. fijl. 



A. The coherent-state propagator 

In Ref. [l?} it is shown that the semiclassical formula 
of the coherent-state propagator, 



K{z 2 ,z l ,T) = {z 2 \e- itlT ' h \z l ), 



(1) 



depends only on complex trajectories of an auxiliary 
classical system governed by the Hamiltonian function 
H(v,u), which is to be built according to the prescrip- 
tion 



H(v,u) 



(z\H\z) 



(2) 



That is, to find H(v,u), one evaluates (z\H\z) and re- 
place z and z* by u and v, respectively. The usual clas- 
sical variables q and p are related to the variables u and 
v through 



1 fq ip 

v/2 U c 



and v 



1 



V2 \b 



ip 



(3) 



where b and c, satisfying be = ft, are related to the vari- 
ances of the coherent state along the position and mo- 
mentum axes. Hamilton's equations written in terms of 
u and v become 



i dH 
ft dv 



and 



i dH 
ft du 



(4) 



Trajectories contributing to the semiclassical propagator 
must satisfy the boundary conditions 



w(0) = z\ and v(T) = z^. 



(5) 



A careful inspection of the dynamical structure defined 
by Eqs. ©-([5]) reveals why the classical variables q and 
p must be complex. Since the boundary conditions given 
by Eq. (|S|) and the evolution time T are both fixed from 
the outset, it is not possible to find in general a classi- 
cal trajectory satisfying that many conditions simulta- 
neously, unless q and p are allowed to be complex num- 
bers. This is the motivation for the change of variables 
(z*,z) -> (v,u). 

Having found the proper trajectory we can evaluate its 
complex action, 



S{zl,z x ,T) 



ih 



(iiv — uv) — H(v, u) 



dt-A, (6) 



where A = f [u(0)v(0) + u(T)v(T)}, and the function 



1 r 

G(Z2,Z!,T) = 2 J o 



1 f 1 fd 2 H{v,u) 
dudv 



dt. (7) 



The semiclassical propagator is then given by 



tiaj . 



e^ 5+e >, (8) 
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where Af = cxp(— ^\z2\ 2 — ^\z\\ 2 ) . Some comments about 
Eq. ((SJ) are in order. First, it is worth mentioning that 
K, is obtained through a quadratic approximation around 
critical paths — the complex classical trajectories — of 
K, expressed in the path integral formalism. Second, it 
is explicitly indicated that, in principle, one should sum 
contributions of all trajectories satisfying the boundary 
conditions. Third, the label z 2 of K is written as z% 
in K- as the trajectories depend only on the value of 
instead of Z2- In the right hand side of Eq. ((HJl, the only 
dependence on Z2 lies in Af '. 

The difficulties to get a semiclassical expression for 
fC* directly from Eq. ((SJ can be better appreciated at 
this point. Contrary to the Van-Vleck propagator, the 
functions S and Q and the classical variables u and v 
are all complex. Then, taking the complex conjugate of 
Eq. ((8j) implies working with the complex conjugate of 
these functions, which although well defined mathemati- 
cally may not offer a straightforward interpretation from 
the point of view of the quantum-classical connection. 

Next we address this issue preserving the mathemati- 
cal structure that was carefully derived and extensively 
discussed in Ref. 1171. 



B. Unified semiclassical formula 



Let us consider now the generic propagator 



K s (z2,Z!,T) = (z 2 \e 



-i£HT/h\ 



(9) 



where H is a time- independent Hamiltonian, £ = ±1, 
and the kets \z\) and \z%) are coherent states. Clearly, 
K-(z2,Zi,T) = Kl(z 1 ,z 2 , T). In this sense, Eq. © con- 
templates propagators and their complex conjugates in 
a unified formula. In addition, we see that K^(z2, zi,T) 
can be obtained from K (z%, z\, T) by means of the change 
H — >• £H [29|]. Furthermore, since £ is nothing but a real 
constant, the mathematical structure previously delin- 
eated readily applies, provided that we consistently em- 
ploy the mentioned change. 

We start our program of implementing the change 
H — > £H with Hamilton's equations (|4]). We get 



u = — — 



h dv 



and 



i£ dH 
h du 



This changes the interpretation of u and v, making their 
roles swap in the dynamics depending on the value of £. 
In order to avoid this issue we define the generalized time 



t e =£t+(l-0T/2. 



(11) 



Explicitly, we see that t+ = t. but for ( = -1 we get 
t- — T — t. This strategy allows us to preserve the 
equations of motion in the same form as Eqs. @, 



in 



i dH 



and 



i dH 



(12) 



where and arc defined by 

i^fo) = u(f(i c )) and w c (i c ) = «(«(*£)), (13) 
with tfa) given by the inverse of Eq. (jllj) . and 



zHz) 



(14) 



In terms of the new functions, the boundary conditions 
given by Eq. ([5]) read 



«e \fi - Of) = 2 i and v «((i + 

or, equivalently, 

u + (0) = z u v+{T) = zl 

U-(T) = zi, V-(Q) = z| 



(15) 



(16) 



We focus now on the functions S^{z^Zi,T) and 
G^(z2 , zi, T), the extended forms of Eqs. © and (J7J. 
Relations (fTB")) give us the rule to rewrite 5^ and in 
terms of u^(t^) and )• We then change the variable 
of integration from t to t^ and the limits of integration 
to (1 - i)T/2 and (1 + £)T/2. Finally, we replace the 
dummy variable t% by t and use the identity 



J F(t) dt = £ j F(t) dt, 



which holds for any F(t) as far as £ = ±1. This procedure 
allows us to write 



d 2 H(v^ H ) 



dt 



(17) 



and 



ih 



dt-A(, (18) 



where A ? = ^ 



For the sake of compact- 



ness of the notation we have introduced the following 
(double) primed variables: 





<4 = 


4 


(10) 




4 



vs(T). 



(19) 



The notation is such that prime (double prime) refers al- 
ways to the initial (final) instant. Notice by (|T9l and (fTS"]) 
that, while u' + = u"_ = z\ and v'! — v'_ = z\, the vari- 
ables u'l , u'_ , v' + , and v"_ arc not fixed by the boundary 
conditions (fT5|) . They are obtained once the solution for 
the trajectory has been found. 

The semiclassical propagator in the coherent-state rep- 
resentation is then finally written as 



traj. 



d 2 Ss 

h dzXdz 1 



1/2 



£(S 5 +6e) 



, (20) 
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where Af = exp(— i|22| 2 — 2"l z i| 2 ) remains unchanged. It 
is worth noticing that for £ = +1 the original formalism 
is fully reproduced. 

Finally, concerning the complex action z\, T) it 

satisfies the following useful relations: 



i 8S4 



i dS-i 



h dz 



i dS- 



Hdv'l ' 



i dS- 



h dz 



hdv' 



i dS + 


i dS + 


h dz\ 


h du' + ' 




(21) 


i dS- 


i dS+ 




h du"_ ' 



and 



r) 1 



(22) 



In addtion, using Eqs. (|2"TT) . the prcfactor of JC^z^, Zi,T) 
can be written as a function of the elements of the tangent 
matrix defined by 



Su'l 
Sv'^ 



= M. 



6v' 



uu w 



One can show that 
i d 2 S + 1 



hdzXdzi My V 



and 



d 2 S^ 



Sv' 



1 



H dzXdz\ Mu 



• (23) 



(24) 



and 



M+ 

UU 



M+ = 



d 2 S 4 



du'dv" dv 



d 2 S+ ( d 2 S+ 



dv'i 



dv'+du'^ 



+ 
-l 



d 2 S+ 
dv'ldu'i 



d 2 s 4 

du''' 



(25) 



M+ = - 



( 9 2 S+ 



2 <? ^'9 2 5 4 



\dv'^_du' + 



du' 



Another set of three equations relating derivatives of S- 
with elements of M_ can be obtained by simultaneously 
replacing +, u and v in Eqs. (|25[) with — , v and u, re- 
spectively. The reason to write the prefactor in terms of 
elements of the tangent matrix is the ease of handling 
them in several situations, specially in numerical treat- 
ments. 

The classical structure we have proposed is such that 
a trajectory (u^(t^) , v^(t^)) is solution of the equations 
of motion in terms of a proper time scale t^. The inter- 
pretation of a forward time evolution from to T is pre- 
served, but while /C+(z|,zi,T) depends on a trajectory 
that propagates from z\ to z|, K-{z 2 , z\,T) depends on 
one propagating from z| to z\. In this sense, comparing 
with the case in which £ = +1, trajectories for £ = — 1 can 
also be interpreted in terms of a backward time evolution, 
which is compatible with the intuition one may construct 
from the exact relation K_ (z 2 , Z\, T) = A' + (z 2 , Z\, —T). 

The set of equations given in this section defines the 
general recipe to obtain the semiclassical version JC^ of 



the exact propagator given by Eq. ((9j). As such, this 
unified formalism constitutes the first important contri- 
bution of this work. All the formal details involved in the 
derivation of original formulas, specially those associated 
to the stationary phase method, can be found in Ref. [TtJ 
for the case in which £ = +1. 



A simple example: harmonic oscillator 

In order to clarify the notation and illustrate the ad- 
equacy of the formalism, we calculate the semiclassical 
version of (z 2 , z\,T) for the harmonic oscillator Hamil- 
tonian, Hh Q . According to Eq. (|T4")) the classical Hamil- 
tonian results 

1 



i?ho(v£,^) = fllO ( + - 



where we have adopted as the coherent-state basis ex- 
actly that one associated to Hho- Then, from Eq. (fl"2"|) 
one gets = iuiv^ and ii^ = —iuju^ whose solutions read 

w c (t c ) = C|e ia, *« and = C|e- iwt «. 

From Eqs. (jT5j) and (|TT|) we get 

Cj = z*e~^ 1+ ^ T ' 2 and C« - z^ 1 ^ 2 , 
so that 

^(t £ ) = z * e *«[t t -(i+c)r/2] j 

u c (t f ) = Zie -MH-{i-t)T/2]_ 
Using these solutions we directly obtain 

G e = — — , S s = ih Zl z 2 e 4 , 

and A^ = ihu^v^ = i/iziz 2 e _IL ^ T '■ The prefactor becomes 

hdz\dz\ i 
The final result is 

fC^(z2,Zi,T) = e~^e~~5 } ~" "~*> e~ 
which is identical to the exact one, 

X 5 (z 2 ,z 1 ,T) = (z 2 |e-^ T ( ata+ 5)|z 1 ) 
= e-^(z 2 \e- i ^ T z 1 ). 



III. SEMICLASSICAL MEASURE OF 
ENTANGLEMENT FOR PURE BIPARTITE 
SYSTEMS 

We now focus on the main task of this paper, namely 
the derivation of a semiclassical measure of entanglement 
for pure bipartite systems via coherent-state propagators. 
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In order to do so we need to extend the results of the 
previous section to bipartite systems. The procedure is 
well known [3(| HH and its generalization for the complex 
conjugate of the propagator is straightforward. 

We consider the coherent-state basis given by |z) = 

\ z xi z y) = \ z x) ® \ z y 

), and the classical variables u ? = 
(u x ,u y ) and = (v x ,v y ). While has a straightfor- 
ward extension, the function requires the change 



dvtduc 



<9 2 iJ(v c ,u c ) d 2 H(v^u 6 ) 



dvidu% 



dvydul 



As far as the prefactor is concerned we need to replace 
the function [i/K){d 2 S^/ dz^dzi) by 



det 



d 2 S e 



' dz% dz\ x dzZ dz\. 



\ dz% dzxx dz. 



d 2 S^ 



det I ^ 



Their scmiclassical approximations are then given by 

(z 1 |p(r)|z 2 ) somi = /C + (zt,z ,T) /C_(z*,z 2 ,T), (30) 

which can then be evaluated by means of complex clas- 
sical trajectories (u^(i^), v^(^)) with specific boundary 
conditions. While, for /C+(zJ, zo, T), the boundary con- 
ditions arc u' + = zo and v" = z*, for /C_ (zjj, z 2 , T), 
they are u" = z 2 and v'_ = Zq. Matrix elements of 
p(T), therefore, can be semiclassically written as func- 
tions of pairs of (generally complex) classical trajecto- 
ries [(u + ,v + ) and (u_,v_)], connected by the fact that 
u' + = z and v'_ = Zq. 

Tracing over the subsystems y, we obtain the matrix 
elements of the reduced density matrix 



[ z lx\Px 



>C + ((z* x ,z* y ),z ,T) 
x lC^z*,(z 2x ,z y ),T)^. 



(31) 



which can be equivalently written as (see Eq. (HH) 



det ' 



(detM+) \ for£ = +l, 
(detM-J" 1 , for£=-l. 



(27) 



Notice that M+ v and M~ u are now 2x2 blocks of the 
tangent matrix. Equations (f2~5j) (and also their versions 
for £ = —1) can also be extended to the case of bipartite 
systems by replacing each second derivative of by a 
2x2 matrix analogously to that of Eq. (|26[) . 

The entanglement of a pure bipartite system composed 
by the subsystems x and y, at the time T, can be quanti- 
fied by the linear entropy of the reduced density matrix, 



SliniPx) = 1 - P{Px)i 



(28) 



where p x = Tr y p and p = \ip(T))(tp(T)\. The purity P of 
the reduced density matrix p x is defined by 



P(px) = Tr x {pl} = Tr x {[Tr yP (T)} 2 } 



(29) 



The information about the entanglement dynamics, en- 
coded in the linear entropy Sn n , is fully contained in the 
purity, which hence is the object of interest in this sec- 
tion. As we are mainly interested on the dynamical be- 
havior of the purity, hereafter we shall denote P(p x ) sim- 
ply by P T - 



A. Semiclassical reduced density matrix 

Assuming an initial state given by |zo) = \zq x ) (8 \zq v ) 
and a generic time-independent Hamiltonian H , the ma- 
trix elements of the density operator in the coherent-state 
representation read 



(zi|p(T)|z 2 ) 



( Zl | e -^/ ft | Zo )(zo|e^/ ft |z 2 ) 
K + (z u z ,T) K_(z ,z 2 ,T). 



To calculate the integral we apply the saddle point 
method [l^llll- The critical points (z*,z y ) satisfy the 
relations 



dz*, 



dz y 



+ tS+{{z* x ,z*),z ,T) 



-\ z v\ 2 + tS-( z oi( z 2x,z v ),T) 



= 0, 



0. 



(32) 



As usual [17( we neglect for it is a low-order term in H. 

According to Eqs. (|2ip . the last equations imply that 
the critical pair of trajectories [(u+,v+) and (u_,v_)] 
contributing to Eq. (|31|) should obey the additional 
boundary conditions u y (T) = z y and v~ (T) = z* . Then, 
given the primary boundary conditions v^(T) = z* and 
u~(T) = z y it follows that among all pairs of trajectories 
contributing to Eq. (|3"Tj) . the critical ones (still complex, 
in general) are those for which the position and momen- 
tum in the y space at the final point are real, having the 
same value for both trajectories, namely, u^(T) — z y and 

Expanding the integrand up to second order around 
the critical pair of trajectories, we get 



{z\x\Px{T)\z 2x ) 



^ M e i\(s++g++s-+g-) 
^ s Vdet M+Vdet M lm 



(33) 

where the bar over the functions indicates that they 
should be calculated with the critical pairs and J\f = 
e -|z | 2 -i|zi*| 2 -i|z 2 *| 2 -|2 B |^ j n addition, 



/ = 



dz*d 



y U*y 



2ni 



exp <^ -5z y Y 5z y > , 



where 



5z , 



([*» 



6 



is the transpose of 5z y , and 



Y 



i d A S- 
-1 



-1 

i d 2 S+ 



The result for the Gaussian integral, 

i d 2 S+ 



I = 



1 - 



hd{z*} 2 

can be alternatively written in terms of the tangent ma- 
trix, 

* d 2 S^ 



h dz 2 



i d 2 S+ 
hd{z* y } 2 



hjM-(M-) hy, 

h£M+ (M+ v )-\, 



where we have defined the column matrix hj,, whose 
transpose reads h^ = (0 1). Using these expressions 
we get 

I 



l-h£M- (M-) \h^M+ (M+) \ 



(34) 



B. Semiclassical Purity 



Now we proceed with the derivation of the semiclassical 
formula for the purity Pt- For convenience, we introduce 
the notation 

(zlx\px(T)\z2x)semi = R (v+ ,u' + ,v'_,u'!_,T), (35) 

where we recall that contributing pairs of trajectories 
(u±(t±),v±(t±)) have boundary conditions = z , 
v'_ = Zq, v" = (z* x ,Zy) and u" = (z 2x ,z y ), and also 
Uy~(T) = z y and v y (T) = z*. The latter conditions state 
that y-position and y-momentum at time T must be real, 
with these two classical quantities defining z y . For the 
sake of clarity, we have eliminated the bar over the tra- 
jectories involved in Eq. (|35[) . 

Noticing that the purity can be written as 

Pt= %x \ Px (T)\w x ), 

we write the semiclassical purity as 
V T = j R(v'l,xi' + y_,n'L,T) 

xR(V"U' V'_ , \J'L , T) d2w * f z * ( 36) 

where the two contributing pairs of trajectories, (u± , v-t) 
and (U±,V±), satisfy, respectively: 

(i) = z ,v'_ = z*,v'L = {w%,zl),u'L = {z x ,z y ), 
u y{T) = z y and v~ (T) = z*; 

(ii) u' + = z ,v'_ =z* ,v; = (z*,w* y ),xs'L = (w x ,w v ), 

U+(T)=w v andK-(T) 



In order to find the critical trajectories [(u±,v±) and 
(U±,V±)] of Eq. (|36|) . we look for its saddle points, 
(w Xl w%) and (z x ,z*). We find the following additional 
conditions: 



V-(T) 
v x (T)-- 



Zn,. i 



ut{T) - 



w x , 

- Zx • 



(37) 





= zo, 


n 


= «>4)> 


< 


v' 


= z 0j 


u" 


= ( Z X , z y)i 


v" 


u' + 


= zo, 


n 


= (K^y)> 


XJ'l 


V'_ 


= zS, 


u'L 


= (W X ,Wy), 





Therefore, all boundary conditions that must be satisfied 
by the critical set of four classical trajectories contribut- 
ing to Vt can be summarized as follows, 



(W X ,Zy), 

( — * — * \ 

(i t) ( 38 ) 

" \<* X : w y ) i 



As discussed previously, the final point of the trajec- 
tory (u + , v+) is connected to the final point of (u_, v_), 
implying the position and momentum in the y direction 
to be real and the same for both trajectories. An equiv- 
alent conclusion applies to (U + ,V + ) and (U_,V_). 
Analogously, in the x direction, we see by Eqs. (|37l) that 
the trajectory (u + ,v + ) is connected to (U_, V_), while 
(u_,v_) is connected to (U+, V + ). This means that the 
trajectories contributing to Vt constitute a set of four 
trajectories whose final and initial conditions are mutu- 
ally connected according to Eq. (|3"5|) . 

A close look at these boundary conditions reveals that 
there exists at least one trivial set of classical trajectories 
satisfying all of them. It corresponds to the trajectory 
starting from u' + = u'_ = U' + = XJ'_ = zq and = 
v'_ = = V'_ = Zq. Hereafter we use tilde to refer 
to this set of four identical trajectories which obviously 
satisfy, in addition, the conditions z x = w x and z y = w y . 

Now a rather important point concerning the con- 
tributing trajectories should be identified. Consider the 
class of time-independent classical Hamiltonians H(v, u) 
deriving from hcrmitian Hamiltonian operators _ff(q, p). 
In these systems a trajectory whose phase-space variables 
are all real at a given instant of time remains real for 
all times [3j|. Since the boundary conditions given by 
Eqs. (|3"5)) assure that the final point is real, the criti- 
cal set of trajectories contributing to Vt has exclusively 
real trajectories. Therefore, once the initial point is com- 
pletely determined there is no other solution to Eqs. (f3"5)) 
but the trivial set discussed above. 

Applying the saddle point method to expand Eq. ([36)) 
around the set of real trajectories we obtain 

V T =lR(v'L, u' + , v'_ , u'L , T) R(V» , U; , V'_ , XJ'L , T) , 

(39) 



where 



with 



d w x d z x 
(2ttz) 2 
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The 4x4 matrix A can be written as 



with 





V c c 



-1 


c c 


+ c b 


c c 


c c 


A a + C a -1 





-1 A b + C b 



where 



A a 

C a 

C b - 

Cc 



hd{u»y 



A b = 



i d 2 S+ 
ft <9K} 2 ' 



i d 2 S- \ I i d 2 S+ 

hdu»du^J ynd{v%y 

i d 2 S+ \ 2 f i d 2 S- 



Kdv»dv») \hd{u''}i 
' i d 2 S+ \ ( i d 2 S- 



h dv"dv'y I \ h du'ldu'y 



D. 



D. 



D, 



and 1/D = 1 - (*j£fp) (|^)- The Ganssian 
integral X then results 

(40) 



I = {[1 - (A, + C a )(A b + C b )]' 

- 2C 2 [1 + (A a + C a )(A b + C b )] + Ct) 



Since the four trajectories are identical, we define 
M uu = M± = M±y , M uv = M± = M± v , 



M vv = M± = M± 



(41) 



x vu — lv - L VU' ^vj. vv — j.vj. vv 

Then, for r or s assuming x or y, 

= h ^ vu M-h s , 



vv 



i d 2 S 4 



(42) 



h^ MuvM" 1 h s , 



with h^ = (1 0). In addition, i?(v" , u' + , v'_ , u" , T) 



R{\" , , V'_ , V"_ ,T)=R, with 



R = 



det M v 



det M uu 

1 - h^M vu (M uu ) \ y h^M uv (M vv ) 



(43) 

Inserting the last results in Eq. (|39|). we obtain 

V T = £~ 1/2 det M uu det M vv (44) 

where 



£ = £ ' 



£" + det M uu det M vv - det A det B 



x det M uu det M vv - det C det D ) - £" 



£' = -4 {det M uu det M vv det A' det B') 
£" = (dctA'VdetBdetD- fdetA'detB'^ 2 



+ 



(det B' 



det A det C, 



A D 
C B 



A' D' 
C' B' 



10 
1 
10 
10 

10 
10 
10 
1 



M, 



(45) 



M. 



Equation fj44[) defines the general recipe for the calcula- 
tion of the semiclassical purity and constitutes, therefore, 
the second important contribution of this paper. A cru- 
cial information emerges from this result, namely, that 
the semiclassical purity strongly depends on the determi- 
nant of sub-blocks of the tangent matrix. This implies the 
purity to be essentially determined by the stability of the 
(real) classical trajectories underlying the corresponding 
classical system. In other words, the semiclassical purity 
is sensitive to whether the trajectory is chaotic or regular. 

Note that Eq. (|4"4"f results 1 for the case of noninteract- 
ing subsystems, in agreement with the result predicted by 
quantum theory. In this case the elements M UrUs , M UrVs , 



and M VtVs , where both r and s may assume x and 



y, with r ^ s, vanish because the subspaces do not cou- 
ple. Then a straightforward manipulation of Eq. (|4"4")l 
leads to the expected result. 

Therefore, given the classical Hamiltonian H (v, u) and 
the center Zq of the initial state, the calculation of the pu- 
rity with Eq. (|44f becomes a problem of classical mechan- 
ics. One may wonder whether the semiclassical formula 
is able to describe the dependence of the purity on the 
characteristics of the initial state other than its centroid. 
However, by examining Eq. (|14[) we realize that iJ(v,u) 
itself has information not only about the physical inter- 
action but also contains quantities that characterize \zq), 
namely its variances b XjV and c x ^ y . 



IV. CASE STUDY: NONLINEARLY COUPLED 
OSCILLATORS 

As an example of application of the formalism we show 
now that, using Eq. (|44|) . the short-time behavior of the 
purity is suitably reproduced. 

Consider the following Hamiltonian 



H = H x (g) l y + l x (g) H y + XH X 



(46) 



where 



H r = 



m r uj 2 q 2 



2m r 



for r — x or y. The initial state |-0o) = \ z 0x) <H> \zo y ) is 
chosen such that \zo r ) is the coherent state associated to 
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H r . The annihilation operator a r and its eigenvalue z$ r 
are 



1 / q r ip r 

a r = — = 1 and z 0r 

y/2 \b r c r 



_L ( _|_ i l^L 

y/2 V b r Cr 



where b r = y/Hj (m r uj r ) and c r = \/m r hiu r . (qa r ,Por) 
gives the location of the center of the wave packet in 
phase space. In terms of the annihilation and creation 
operators the Hamiltonian is written 

t; 



H 



ffllyOL^OLy 



hTa x a x a v a y + e , (47) 



Kl x a' x a x 

where f2 r = u r + T/2, T = Xhoj x ujy, and eo = h(uj x + 
uj y )/2. According to (|14p the underlying classical Hamil- 
tonian is 

iJ(V, U) = htt X V X U X + MlyVyUy + HTV X U X VyUy + Eg . 

The classical trajectories can be readily integrated and 
are given by 



/ u x (t) 

Uy(i) 

v x (t) 
V Vy(t) 



e -Kt ^ 



u y e 
vLe +Kt 



(48) 



where = i (Tl x + Tu' v' ) and X y = i (tt y + Yu' x v' x ). 
The tangent matrix can be written as the product of two 
matrices Mi and M 2 such that 



(49) 







Su' x \ 


5u>> 




H 


K 


= M 2 Mi 




V K J 







where 



Mi 



/ 


1 











-au' y v' x 


1 












av' x v'y 


1 




\ 







au£i4 


1 



with a = «TT, and 
( 



M 2 = 



V 



-A X T 









3 — XyT 

















= +A v T 



As we are interested just in the trajectory starting from 
u' = Zo and v' = Zq, we obtain 

dctM uu = e-^+^ T (l-a 2 \z 0x \ 2 \zo y \ 2 ), 
detM vv = e+^+ A *> T (l-a 2 |z 0:c | 2 K| 2 ), 

and, according to Eq. (|4"5| . 



dctA = a 2 |zox| 2 (^) 2 e- (A -- A « )T , 

dctB = a 2 \z 0x \\z 0y Ye+^~^ T , 

detC = a 2 (z *J 2 |z 0y | 2 e+(^- A ^ T , 

dctD = a 2 {z 0x ) 2 \zo v \ 2 e^^- x y^ T , 

detA' = azQ X ZQ y , 

detB' = — azoxZoy. 



Then, Eq. (J44J) becomes 



v /l + 6x + a; 2 (3 + 2a;) 2 ' 



(50) 



where x = IzoaiPkoj/pr 2 ! 12 . It is important to compare 
the semiclassical result with the exact one. Using com- 
mon techniques of the quantum formalism we obtain 



p. L, = e -2|»„.| : 



E 



\ZQx 



2(n+m) 



-4\z 



2 rrr(n-m)i 



J 1 7T) I 



(51) 

This result is clearly different from the semiclassical one. 
In particular, we see that, at the instant 27r/r, the quan- 
tum result predicts the total recoherence of the subsys- 
tem, i.e., P27r/r = 1- The semiclassical formula, however, 
results in a monotonically decreasing function of T, and 
as such is not able to reproduce the recoherence. On the 
other hand, let us consider the short-time scale (FT < 1). 
In this case the quantum result simplifies to 



Pi 



1 



\z Qx \ 2 \z 0v \ 2 T 2 T 2 



(52) 



which accurately agrees with the semiclassical purity Vt 
in this limit. 

This example revealed the limitations of our semiclassi- 
cal formula. We see that the approach is not able to cap- 
ture the physics of recoherences, which is associated with 
important quantum phenomena such as interferences and 
revivals. Actually, we have seen that the semiclassical pu- 
rity reproduces accurately the exact result only within a 
very short time scale (TT < 1) that is much shorter than 
the one in which recoherence occurs (TT = 2%). It has 
been shown in the literature that interference phenomena 
can be reproduced by semiclassical approaches involving 
more than one trajectory (see, for instance, Ref. [34]]). We 
expect a similar strategy to be able to improve our results 
for longer times. This, however, requires one to carefully 
revisit the formalism looking for further contributing tra- 
jectories, which do not obey exactly the boundary con- 
ditions ([55)1 . 

As far as the entanglement is concerned — here mea- 
sured by the linear entropy Sn n = 1 — Pt — we may 
write 



2 Hj n t Tn- T„ 



(53) 



where we have defined the dimensionless time T r = ui r T . 
In this expression Hi nt corresponds precisely to the clas- 
sical version of the interaction Hamiltonian given in 
Eq. (|4"r?)) . Notice that the short-time entanglement grows 
proportionally to the magnitude of the interaction, as 
expected. Surprisingly, however, it does not depend on 
h at all, thus corroborating our claim that the onset of 
the entanglement dynamics can be described in terms of 
classical mechanisms. 

It is worth noticing that our semiclassical formula does 
predict a dependence on h, in general. Consider, for 
instance, an arbitrary classical function H(q 7 p). The 
application of usual quantization rules to this function 
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(see, e.g., Ref. J3j|) produces an h- independent opera- 
tor H(q,p). However, the classical Hamiltonian enter- 
ing in our recipe is given by (z\H(q,p)\z) — H(q,p) + 
En>o^"/»(9'ri' which generally depends on h [17|- It 
follows that the stability matrix and the semiclassical pu- 
rity will depend on h as well. However, in the regime of 
large actions and energies this dependence manifests as 
a perturbation to the dynamics generated by H(q,p) so 
that our claim remains valid. 



V. FINAL REMARKS 

We have derived a semiclassical formula for the purity 
of pure bipartite systems initially prepared in a prod- 
uct of coherent states. Since we are here concerned only 
with pure states, our formula turns out to be a direct 
semiclassical measure of entanglement. As a preliminary 
step towards the development of our formalism, we have 
derived a unified semiclassical formula which is able to 
approach both propagators and their complex conjugate. 

Our result for the semiclassical purity is given in terms 
of a very compact formula (j4"4")) which is shown to de- 
pend only on the trajectories of an auxiliary classical sys- 
tem. Specifically, the short-time entanglement dynamics 
is proven to depend exclusively on the elements of the 
tangent matrix, which defines the local stability of the 
classical trajectories. As a consequence the initially sep- 
arable wave functions get spread and then entangle ac- 
cording to a rate that strongly depends on whether the 
corresponding classical trajectory is chaotic or regular. 

Finally, in order to illustrate the theory, the formal- 
ism has been applied to the problem of two nonlincarly 
coupled oscillators, whose dynamics is rich in quantum 



effects, such as collapses and revivals. The semiclassical 
approximation has shown to reproduce exactly the entan- 
glement dynamics in the short-time regime. This is con- 
sistent with the approximations underlying the method. 

Our results are in qualitative consonance with those 
reported in Refs. QJ, HU and give additional analyti- 
cal support to the widely known fact that the entangle- 
ment dynamics in the regime of short times depends on 
the characteristics of the classical point in phase space 
in which the initial state has been centered (see, e.g., 
Ref. @). Moreover, they emphasize the fact that the 
short-time entanglement is promoted essentially by clas- 
sical mechanisms, which here have been identified to be 
the stability of underlying classical trajectories. Our find- 
ings provide, therefore, analytical support to the numer- 
ical results of Refs. @, H| which show that it is possible 
to mimic the entanglement dynamics in terms of cntropic 
measures defined in the Liouvillian theory. 

The natural continuation of this work consists in 
extending the formalism to spin degrees of freedom. 
Moreover, even though we have assumed the initial state 
to be a product of coherent states the generalization 
of the semiclassical purity to arbitrary initial states is 
possible. Research on these topics are now in progress. 
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